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Abstract 

We calculate the diffuse high energy (TeV - PeV) neutrino emission from hyperflares 
of Soft-Gamma Repeaters (SGRs), like the hyperflare risen from SGR 1806-20 on 
December 27 of 2004, within the framework of the fireball model. The fireball model 
for gamma-ray bursts (GRBs) can explain well the main features of this hyperflare 
and the subsequent multi-frequency afterglow emission. The expected rate, ~ 20 — 
100 Gpc _3 day _1 , of such hyperflares is well in excess of the GRBs rate. Our result 
shows that the contribution to the diffuse TeV-PeV neutrino background from such 
hyperflares is less than 10% of the contribution from GRBs. We also discuss the 
high energy cosmic rays (CRs) from these sources. 
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1 Introduction 

It is a general consensus that relativistic shocks can accelerate nuclei to very 
high energies through the Fermi acceleration mechanism. By interacting with 
photons or baryons, such high energy nuclei (mostly protons) can generate 
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pions and the latter cascade into neutrinos and leptons. So, for a relativistic 
jet that sweeps up its surrounding medium, it could be an efficient high-energy 
neutrino producer, if the jet is also surrounded by intense photon fields. Three 
famous kinds of such astronomical objects are gamma-ray bursts (GRBs; Wax- 
man 1995; Vietri 1995; Waxman & Bahcall 1997; Bahcall & Meszaros2000; 
Meszaros& Waxman 2001; Guetta, Spada & Waxman 2001; Dermer & Atoyan 
2003; Asano & Meszaros2008), active galactic nuclei (Axford 1994; Atoyan & 
Dermer 2001; Dermer, Ramirez- Ruiz & Le 2007; Berezhko 2008; Abbasi et al. 
2008), and micro-quasars (Levinson & Waxman 2001; Distefano et al. 2002), 
all of which are usually discussed as high-energy neutrino and cosmic ray (CR) 
sources. For some thorough reviews on astrophysical neutrinos and their con- 
nection to CRs, we would like to refer the readers to the references Halzen 
et al. (2002) and Becker (2008). In this paper, we would focus on another 
kind of such astronomical objects, i.e., soft-gamma repeaters (SGRs) which 
are widely accepted as magnetars (pulsars with super strong magnetic field 
~ 10 15 G; Duncan & Thompson 1992). Zhang et al. (2003) proposed a model 
for neutrino production by magnetars in their steady phase of periodic emis- 
sion. In contrast, some short neutrino bursts could also be produced by a 
relativistic outflow during the violent giant-flare phase of SGRs. 

Giant flares are distinguished from common SGR bursts by their extreme en- 
ergies (~ 10 44 ergs) emitted during their initial short (~ 0.1 s — 0.5 s) pulses 
of soft gamma rays followed by subsequent emission lasting hundreds of sec- 
onds showing pulsations associated with the spinning neutron star (Woods & 
Thompson 2006). Among the observed three giant flares during the last four 
decades, the brightest one originated from the SGR1806-20 on 2004 December 
27 (Hurley et al. 2005; Palmer et al. 2005; Terasawa et al. 2005; Mazets et al. 
2005) has an energy release exceeding 10 46 ergs, which is two orders of mag- 
nitude higher than the energy release of the other two. The huge difference in 
the luminosity urges us to consider this kind of events separately and name 
it "hyperflare" (Popov & Stern 2006). The mechanism triggering hyperflares 
remains a matter of debate. Lugones (2007) proposed a model in which the 
core conversion of an isolated neutron star with a magnetic field of ~ 10 12 G 
and a fallback disk around it can trigger a giant flare. In the popular mag- 
netar model, giant flares result from a global magnetic rearrangement within 
the crust of the magnetar (Thompson & Duncan 2001). Since giant flares have 
many similarities to GRBs which can be well understood by the fireball model 
(Piran 1998; Meszaros2002), different scenarios have been proposed within the 
frame work of the fireball model (Wang et al. 2005; Dai et al. 2005; Yamazaki 
et al. 2005), to explain the abundant multi- frequency afterglow data of this 
hyperflare (Gaensler et al. 2005; Cameron et al. 2005; Taylor et al. 2005) and 
the flare itself might be the emission from the internal shock and/or the photo- 
sphere of the fireball if the relativistic outflow is variable (Nakar, Piran & Sari 
2005; Ioka et al. 2005). Ioka et al. (2005) estimated the number of high-energy 
neutrino events from this hyperflare and argued that the neutrino flux should 
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be detected by current neutrino observatories such as AMANDA (Ahrens et 
al. 2002; Halzen, Landsman & Montaruli 2005), which may put constraints on 
the flare mechanism. Fan, Zhang & Wei (2005) considered the production of 
neutrinos with typical energy of 10 14 eV through photomeson interaction of 
X-ray tail photons with ~ 10 17 eV CRs accelerated in the external forward 
shock by this hyperflare. They found that the neutrino fluence produced in 
the external shock is too weak to be detected. 

The rate of hyperflares is very uncertain because of low statistics based on 
a small sample. Their intrinsic low release energies relative to GRBs make 
them impossible to be detected beyond ~ 30 — 40 Mpc by BATSE and ~ 
70 Mpc by Swift (Hurley et al. 2005). Lazzati et al. (2005) gave a slightly 
less stringent limit for the rate of hyperflares: < 1/130 yr" 1 per a Milky- 
Way- like galaxy. Popov & Stern (2006) argued that the rate is ~ 10~ 3 yr _1 . 
They further conservatively estimated the expected rate of hyperflares to be 
~ 20 — 100 Gpc~ 3 day -1 (Popov & Postnov 2007). Recently, Lorimer et al. 
(2007) reported a discovery of a strong millisecond extragalactic radio burst 
with peculiar properties and estimated that the cosmological rate for this radio 
burst is ~ 50 Gpc~ 3 day -1 which is in correspondence with the statistically 
estimated rate of hyperflares. Meanwhile, the millisecond time scale of the 
radio burst is consistent with an event in the magnetosphere of a magnetar, 
indicating that both the millisecond extragalactic radio burst and hyperflares 
may possibly come from the same source: magnetars (Popov & Postnov 2007). 

The above rate of hyperflares in the universe is well in excess of the GRBs 
rate, while the total energy release of a hyperflare is much lower than that of a 
typical GRB. Based on the estimated rate and the fireball model proposed by 
Ioka et al. (2005) for the hyperflare of SGR 1806-20 (see Section 2) in which 
relativistic protons are accelerated by an internal shock and target photons 
are dominated by the hyperflare, we describe the neutrino production process 
within hyperflares in Section 3. We estimate the diffuse TeV-PeV neutrino flux 
from hyperflares, compare our results with that from magnetar steady phase 
and from GRBs in Section 4. We also give a discussion on high-energy CRs 
from hyperflares in this section. We summarize our results and conclusions in 
Section 5. 



2 The fireball model for the hyperflare from SGR 1806-20 



The hyperflare of SGR 1806-20 is the only event up to now and its spectrum 
may be either thermal (Hurley et al. 2005) or non-thermal (Mazets et al. 
2005; Palmer et al. 2005). The baryon load in the outflow is less constrained 
in contrast to GRBs, then the outflow produced the hyperflare may be either 
baryon-poor or baryon-rich (Ioka et al. 2005). Based on the fireball model, we 
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describe these two scenarios as follows. 

The isotropic soft 7— ray energy of the December 27 event, E 1 ~ 3 x 10 46 -E 7)46 .5 
ergs, released within a time interval of to ~ 0.1i ,-i s from somewhere near the 
surface of the magnetar with radius of r ~ 10 6 cm, could create an optically 
thick fireball with an initial temperature of 

To~(-^) 1/4 ~ 300 ^ 4 4 , 5 r -: 6 1 / 2 keV, (1) 

where a is the radiation density constant, L ~ £ 7 /£ 7 ~ 3 x 10 47 L , 4 7.5 ergs 
s _1 and £ 7 is the conversion efficiency of total energy into gamma-rays (Ioka 
et al. 2005). The subsequent fireball evolution depends on the dimensionless 
entropy rj — L / Mc 2 which has a critical value (Meszaros& Rees 2000) 

/ L a T \i/4 1/4 _ 1/4 

^ = (i^r) ~ 100 ^47, -0,6 • (2) 

If 77 < 77*, the fireball is baryon-rich, then the photosphere appears in the 
coasting phase and almost all the energy goes into the kinetic luminosity 
of the outflow, L kin ~ L . While if 77 > 77*, the fireball is baryon-poor, the 
photosphere appears in the acceleration phase and then a small fraction of 
the energy goes into the kinetic luminosity of the outflow L k m ~ L 7 x 77^/77 
(Ioka et al. 2005). The above two scenarios are both possible for the December 
27 hyperflare. The particular parameters adopted for the two scenarios by Ioka 
et al. (2005) are listed below 

Baryon Poor: 77 ~ 10 4 , L kin ~ 10- 2 L 7 ~ 10 45 - 5 ergs s _1 , V ~ 100, At ~ 10" 4 
s, 

Baryon Rich: 77 ~ 10, L kin ~ 10L 7 ~ 10 48 - 5 ergs s" 1 , T ~ 10, At ~ 10" 1 s. 

where V is the Lorentz factor of the outflow and At is the variability timescale 
of the hyperflare. It should be noted that the photosphere emission is ther- 
mal. If the internal shocks occur above the photosphere as in the baryon-rich 
model, the non-thermal shock emission will dominate the photosphere ther- 
mal emission if 77 < lOO?]*^^ 8 , where ^ s is the conversion efficiency of kinetic 
energy into photons (Ioka et al. 2005). 



3 The neutrino production in hyperflares 

The photon distribution in the comoving frame of the internal shocked region 
is usually assumed to be isotropic and described by monoenergetic number 
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density dn 1 /de 1 at energy e 7 . The protons are usually assumed to be acceler- 
ated to a power law distribution by internal shocks, dn p / de p oc e p 2 where the 
proportional coefficient l/ln(e Pimax /e Pimin ) as the fraction of the total energy 
that is contributed by each decade of energy is about 0.1. When bathed in 
the hyperflare photon field, the high-energy protons would lose their energy 
by jry interaction, resulting in plenty of pions. Due to the pion production, 
the fractional energy loss rate of a proton with energy e p = •y p m p c 2 is (Stecker 
1968; Waxman & Bahcall 1997) 




where <7 P7 (e 7 ) is the cross section of the photopion interaction for a target 
photon with energy e 7 in the proton rest frame, k P7 is the inelasticity coefficient 
defined as the average fraction of energy lost to the pion and e 7)t h = 0.15 GeV 
is the threshold energy. Using the A-resonance approximation, Eq. (3) can be 
casted into 

V, °> j W (4) 

e A /27p 



where <7 pea k — 5 x 10~ 28 cm 2 and K pcak ~ 0.2 are the values of a and R at 
e 7 = ca = 0.3 GeV where the cross section peaks due to the A resonance, and 
Ae ~ 0.2 GeV is the peak width. 

Under the A-resonance approximation, the neutrino spectrum is totally de- 
termined by the parameter / P7 ~ r s h/rct P7 , the fraction of energy lost by 
protons to pions, where the internal shock radius r s h is ~ 2r 2 cAt. At each 
interaction, a proton loses ~ 20% of its energy which is distributed roughly 
equally among the products of the decay processes 7r ± — > /i ± + v^ipp) — > 
e ± + z/ e (z/ e ) + + v^. In addition the neutrino oscillation will change neutrino 
flavor from v e : : v T = 1 : 2 : at the source site to 1 : 1 : 1 at the earth, 
which will reduce the observed muon neutrino flux by a factor of 2. Once we 
know the spectrum of the target photon, we can calculate the monoenergetic 
muon neutrino flux as 

$,, P7 = min(l, /p7 )i i ^-, (5) 



where d is the luminosity distance of the source, ^ is the fraction of the 
total kinetic energy converted into accelerated protons and e Uttyp is the typical 
energy of neutrinos produced by the p7 interactions. Below we consider two 
different types of photon spectra of SGR hyperflares. 



5 



3. 1 The thermal spectrum of hyperflares 



Hurley et al. (2005) reported that the observed energy spectrum of the hard 
spike of the December 27 hyperflare is consistent with a cooling blackbody 
with average temperature ksTo^ = 175 ± 25 keV, and thus the photon peak 
energy in the observer frame is e° b p Cak ~ 2.7A;^T b s - The differential photon 
field density in the comoving frame of the internal shocked region is therefore 
dn 1 /de 1 = (8ir/h 3 c 3 )(e 2 /(e e '< /kBT - 1)), where T = T obs /r and e 7 = e° bs /r. 
From n 7 = / °° ^de 7 ~ L 7 /(4vrr 2 cre° b p s cak ), we have (Ioka et al. 2005) 



J 7 x 10 18 L 7i47 .5 e 7 5 3 I7 5 AtZl cm" 3 (Baryon Poor) 
n7 ~|7x 10 17 L 7 ^ 47 .5 e~\ 3 TZ 5 Atzj cm" 3 (Baryon Rich) ' [ > 

where e 7 , 5 . 3 = e° b p s cak /200 keV. 

Inserting the differential photon number density into Eq. (4), we obtain the 
analytical form of the fractional energy loss rate as 



_ x C/ 7 c Ae 2 

^P7 ~ Cpeak^peak^ X X 

£7, peak ^A 



ln(l - e -«) , (7) 



where the dimensionless parameter x = 6a/ '(27 p /c_bT) and the photon energy 
density U 1 = n 7 e 7)Pca k = L 7 /(47rr 2 cr 2 ). Then / P7 can be derived analytically, 



f P1 ,th^x ln(l - e 



j 0.28 L 7j47 .5 Cia r 2 - 4 At:} (Baryon Poor) 

\ 2.8 L 7i47 . 5 e-i 3 rr 4 AtZ\ (Baryon Rich) ' ^ 

which peaks at rc pca k = 2.7/tr 2 eA^ P c 2 /(8e° bs e° b p Cak ) 1.8. Approximation 
of the scaling function |ln(l — e~ x )\ yields e~ x for x > :r pca k, and |lnx| for 
j; < 1. Thus, the exponential suppression of the neutrino spectrum appears 
for e° bs < 3/tr 2 eAm p c 2 /(16e° b p Cak ), and high energy region of neutrinos scales 
as x 2 |lnx|. Meanwhile, we get the peak energy of the neutrino spectrum 



obs 3RT 2 e A m p c 2 ( 560 T 2 . eZ\ 3 TeV, (Baryon Poor) 

^peak " ' 16e cbs cak ~ j 5.6 F 2 £ -l ^ TeV) ( Bary on Rich) 
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3.2 The non-thermal spectrum of hyperflares 



The photon spectrum of the 2004 December 27 hyperflare is also likely to 
be non-thermal, making this event resemble a short, hard gamma-ray burst 
(SHB). Palmer et al. (2005) argued that such an energetic SGR flare may 
indeed form a subclass of GRBs. With a luminosity of ~ 10 47 ergs s -1 , such 
flares can be detected by BATSE up to ~ 40 Mpc, suggesting that a consid- 
erable fraction of the BATSE SHB sample is compromised of similar extra- 
galactic hyperflares (Palmer et al. 2005; Hurley et al. 2005; Popov & Stern 
2006; Nakar 2007). Another two short hard GRBs, 051103 and 070201, are 
recently reported to show evidences as hyperflares from SGRs in the nearby 
M81 and M31 galaxies, respectively (Frederiks et al. 2007; Mazets et al. 2008). 
So although very uncertain, at least part of hyperflares from SGRs would have 
a similar spectrum to SHBs. Since the non-thermal component is negligible 
in the baryon-poor model, we only consider a typical non-thermal hyperflare 
spectrum in the baryon-rich scenario with a break energy e° b fe s ~ 200 keV 



(10) 



where e 7) & = e° b b s /T and rib ~ U 1 /2e 1 ^ 2 - After performing the integration in 
Eq. (4), we can approximate the energy loss rate of protons by (Waxman & 
Bahcall 1997; Murase & Nagataki 2006) 



t-'-^a „« ^xHrf^) 1 ^' ei< ^ b ill) 

r p7 ~~ o °pcak«pcak x S / _ N 1-/3 , [11) 

and thus 



obs \ Q — 1 i i 

v \ ! obs ^ obs 

obs J 1 t / h 



/p 7 ,non-th - 1-4 L 7 ,47.5 e 7 ,5.3 r i ^-1 { , "<L \ a-\ u u ( 12 ) 



r° bs V 



^obs ^_ ,obs 
e u > e u,b 



where e° bs ~ RT 2 t\m p c 2 / (8e° b 6 s ) pa 4Tf TeV is the neutrino break energy in 
the observer frame. The power law neutrino spectrum resulting from a non- 
thermal target photon field is quite different to the one predicted from a 
thermal photon field as in Eq. (8). 
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4 The diffuse high-energy neutrinos and CRs from hyperflares 



Similar to GRBs, a high-energy neutrino flash accompanying an SGR hyper- 
flare would be expected if the source is nearby and energetic. Below we discuss 
the probability of detecting TeV-EeV muon neutrinos from the SGR 1806-20 
hyperflare by IceCube, using the following formula (Razzaque et al. 2004; Ioka 
et al. 2005) 

P(e„) = 7 x 10^ 5 (e,/10 4 - 5 GeVf, (13) 



where (3 = 1.35 for e u < 10 4 5 GeV while (3 = 0.55 for e u > 10 45 GeV. The 
number of muon events from muon neutrinos above TeV is given by 

PeV 

iV M = A dct t J P{e v )^ v ^de v , (14) 

TeV 



where the geometrical detector area of IceCube A dct ~ 1 km 2 (A dct ~ 0.03 km 2 
for AMANDA). Then the number of muon events from p7 neutrinos at typical 
energy (e^typ = e^ b p Ca k f° r thermal spectrum and e^ typ = e° b b s for non-thermal 
spectrum) is 



A de lt()P (6;y ) typ)Cl/,typ < ^ ) ^,p7 

!1.3 min(l, / P7 ,th)& -i£kin,45.5dr 2 (e 7 ,5.3rl)^ _1 Adet,i*o ,-i (Baryon Poor) 
2500 min(l,/ P7ith )e i ,-iL kini 48.5rfr 2 (e7,5.3r 2 ) /3 - 1 A dctil to,-i (Baryon Rioflip) 
2200 min(l, fp*f t non~ th)£i,— l-^kin,48.5 ^'(rD^^det,!^,-! (Baryon Rich) 

where the distance of SGR 1806-20 is taken to be 10di kpc (see Ioka et al. 
2005; Dai et al. 2005; however d ~ 15 kpc in Corbel & Eikenberry 2004). 
Due to the model uncertainties, the number of muon events spans more than 
three orders of magnitude. To see clearly the behavior of the number of muon 
events depending on the model parameters, we can fix the values of £j, d and 
set At free in the baryon-poor model, then the event number is approximately 
proportional to 77*/?7(min[?7, 7]^}) 2f3 ~ 2 . Parameter space in the 77 — At plot where 
IceCube can detect more than one muon event from SGR 1860-20 is shown in 
Fig. 1 of Ioka et al. (2005). It can be seen that a baryon-poor outflow with en- 
tropy less than 2500 can trigger one event. While if the outflow is baryon-rich, 
IceCube can detect about one event even when a hyperflare considered here is 
located ~ 200 kpc away. If such detection comes true, it will provide indepen- 
dent evidence for the picture in which relativistic outflows produce hyperflare 
electromagnetic emission and constrain the parameters of the baryon loading, 
the bulk Lorentz factor of the fireball, the efficiency of energy conversion and 
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the variability timescale of the hyperflare. It should be noted that we neglect 
the neutrino production through the pp reaction in this paper. 

Now we turn to the observations. The AMANDA-II detector was running 
to search for down-going muons cascading into high-energy gamma-rays and 
neutrinos, when the hyperflare on December 27, 2004 saturated many satellite 
gamma-ray detectors. However, the data revealed no significant signal which 
put an upper limit on neutrino flux of the hyperflare: < 0.4(6.1) TeV m -2 s~ ] 
for an energy spectrum index —1.47 (—2) (Achterberg et al. 2006). This limit 
would suggest a baryon-poor outflow, which means that the diffuse neutrino 
flux from the hyperflare is much lower than 1.0 x 10" 10 GeV cm -2 s -1 , far 
below the IceCube sensitivity ~ 8.0 x 10~ 9 GeV cm" 2 s" 1 after one year 
(Ahrens et al. 2004; Hoshina et al. 2008; Halzen 2008). But IceCube, cur- 
rently under construction at the South pole, can potentially detect TeV-PeV 
neutrinos from a single hyperflare originating from a baryon-poor outflow (see 
Eq. fl~5]) . Another detector Km3NeT, a planned experiment in the Mediter- 
ranian Sea to complement the IceCube, will have a better detectability for a 
south-hemisphere source like SGR 1806-20 due to its low background of atmo- 
spheric muons (Katz 2006). As estimated above, in the baryon-rich model, a 
hyperflare with the same energy as the hyperflare of SGR 1806-20 can be de- 
tected by IceCube within ~ 200 kpc. If AMANDA-II really had not observed 
high energy neutrinos from SGR 1806-20, then IceCube can only detect hy- 
perflares in our galaxy. SGR 1806-20 and SGR 1900+14 showing giant flares 
are associated with massive star clusters (Fuchs et al. 1999; Wachter et al. 
2008), and SGR 0525-66 is associated with a supernova remnant (Gaensler 
et al. 2001; Eikenberry 2003). The best sites for IceCube to search for high 
energy neutrinos from SGRs are active star-formation regions in our galaxy, 
such as Westerlund 1, a young massive star cluster, in which several neutron 
stars and magnetars were discovered (Muno et al. 2006; Muno et al. 2007; 
Clark et al. 2008). For electromagnetic signals, the detection distance can be 
as far as ~ 70 Mpc by the Swift's Burst Alert Telescope (BAT), which gives 
an excellent opportunity to observe extragalactic giant flares and hyperflares 
from SGRs (Hurley 2005). Popov & Stern (2006) proposed the most promising 
targets for such observations, e.g., Virgo Cluster, NGC3256 etc. Although as a 
smaller detector, the Fermi's Gamma-Ray Burst Monitor (GBM) is less sensi- 
tive than BATSE and thus its detection distance for hyperflares is < 40 Mpc, 
can augment the Swift's high energy sensitivity. Via the Fermi's Large Area 
Telescope (LAT), the detection can be naturally extended to the GeV region 
(Gehrels et al. 2007). The broad-band studies will enable the investigation of 
the hyperflare spectral and temporal evolution over 7 orders of magnitude in 
wavelength, shedding light on the mechanism of hyperflares. The gamma-rays 
decayed from the neutral pions produced by hadronic interactions have a flux 
and energy comparable to neutrinos. These TeV-PeV gamma-rays might be 
detected by the southern Cherenkov telescope, e.g., H.E.S.S, if the gamma-rays 
can escape the pair production absorption (Kohnle et al. 2000). 
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Besides possible detection of neutrinos from a single hyperflare, the contri- 
bution to the neutrino background from such hyperflares is also expected to 
be important since the hyperflare explosion rate is high. The total energy 
spectrum of the accelerated protons can be written as 



Then we can estimate the diffuse neutrino background flux from hyperflares 
based on the template of the SGR 1806-20 hyperflare by 



where H = 71 km s _1 Mpc -1 , /& is the beaming factor, .Rhf(O) is the to- 
tal hyperflare rate at z — 0, and f z is the correction factor for the possible 
contribution from high redshift sources. Although we can always choose At 
to set that min(l,/ p7 ) = 1, the diffuse TeV-PeV neutrino flux, however, also 
depends on the baryon loading of the outflow, the explosion rate, the efficiency 
of energy conversion and the geometry-corrected energy of the hyperflare, and 
therefore it is more uncertain than the high energy neutrino flux from a sin- 
gle source. Among these four parameters, the efficiency £j is more stringently 
constrained because it is less dependent on the model of hyperflares, the other 
three parameters have large uncertainties due to low statistics with using only 
one sample. In addition to strong hyperflares, SGRs also show less intensive 
but more frequent giant flares, e.g., the giant flares of March 5, 1979 from 
SGR 0525-66 with an isotropic energy E iso = 3.6 x 10 44 ergs, and Aug 27, 
1998 from SGR 1900+14 with E iso = 5.2 x 10 43 ergs. Their galactic rate is 
estimated to be about 0.05 — 0.02 yr" 1 (Woods & Thompson 2006), one or- 
der larger than that of hyperflares. The less energetic giant flare may have a 
larger beaming factor than typical hyperflares, thus the giant flares and hy- 
perflares would possibly provide the same contributions to the diffuse neutrino 
background. Due to the low statistics, it is unknown whether giant flares and 
hyperflares from SGRs form a continuous luminosity distribution. We cannot 
estimate the contribution of high energy neutrinos from hyperflares with en- 
ergy larger than the December 27, 2004 event. Nevertheless, the baryon- rich 
model with high explosion rate and large geometry-corrected energy predicts 
high diffuse neutrino flux. To be specific, we estimate a TeV-PeV diffuse neu- 




(16) 



e 2 J. 




x f 0.86 x 10 14 & _ .5 £kin,45.5 *o,-i (BaryonPoor) 
\ 0.86 x lO" 11 & _o. 5 £ k in,48.5 *o,-i (BaryonRich) 
GeV cm" 2 s" 1 str" 1 



(17) 
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trino flux of 4 x 10~ n GeV cm~ 2 s _1 str -1 by using an optimistic set of 
model parameters (rj = 10, £j = 0.5, fi> = 0.25, L 7) 47.5 = 1.0, to — 0.1 and 
i?HF = 100 Gpc~ 3 day -1 ). This value is about three orders of magnitude higher 
than that estimated in magnetar steady phase (Zhang et al. 2003), while more 
than one order of magnitude lower than that from high-luminosity GRBs and 
low-luminosity GRBs (LL-GRBs) (Waxman & Bahcall 1997; Murase & Na- 
gataki 2006). 

In addition to be served as high-energy neutrino factories, hyperflares are 
also possible sources of CRs (Asano et al. 2006). A baryon-rich outflow and a 
baryon-poor outflow can accelerate protons to the maximal energies of 7x 10 18 
eV and 4 x 10 16 eV, respectively (Ioka et al. 2005). Recently, Wang et al. 
(2007) and Budnik et al. (2008) argued that the trans-relativistic supernovae 
accompanied by LL-GRBs may be the mecca of CRs within the energy range 
10 17 — 10 19 eV due to their energetic outflows and high explosion rate. Is it 
possible that hyperflare are important sources of CRs in this energy range? For 
a baryon-poor outflow, low magnetic field can not accelerate CRs to sufficient 
high energies and there is no adequate kinetic energy available for CRs. While 
for a baryon-rich outflow, the kinetic energy is ~ 10 47 ' 5 ergs, four orders of 
magnitude lower than that of LL-GRBs. On the other hand, the explosion 
rate of hyperflares is two orders of magnitude higher than the latter. So even 
if hyperflare have a baryon-rich origin, their contribution to CRs within the 
energy range 10 17 — 10 19 eV is only ~ 1% of that of LL-GRBs. 



5 Conclusions 



In this paper, we have calculated the diffuse TeV-PeV neutrino flux from 
hyperflares of SGRs based on the template of the famous 2004 December 
27 event and the fireball model. Because of a possible high explosion rate, 
albeit very uncertain, hyperflares have a contribution, although likely less 
than 10% of the contribution from GRBs, to the diffuse high-energy neutrino 
background. However, the nondetection of the neutrinos from SGR 1806-20 by 
AMANDA II indicates that the hyperflare may have a baryon-poor origin and 
thus its contribution to neutrino background can be negligible. It is hard for 
IceCube to detect the diffuse high energy neutrinos from hyperflares, though 
which can put severe constraints on the hyperflare explosion rate. However 
the neutrino detectors of IceCube and Km3NeT can optimistically capture 
the TeV-PeV neutrinos from a single event like the hyperflare by SGR 1806- 
20, which may provide independent clues on the trigger mechanism of the 
hyperflare in addition to its electromagnetic emission. More observations are 
required to improve the statistics and verify what scenario is correct. 
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